Multi-Chimera States in FitzHugh-Nagumo Oscillators 
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We demonstrate the existence of chimera states in a ring of identical oscillators described by 
FitzHugh-Nagumo equations with nonlocal coupling. This class of elements serves as a paradig- 
matic model in neuroscience, chemical oscillations, and nonlinear electronic circuits. Applying a 
phase-reduction technique we show that off-diagonal nonlocal coupling is a crucial factor for the 
appearance of chimera states, which consist of coexisting domains of coherent (phase-locked) and 
incoherent oscillators. Surprisingly, we find that for increasing coupling strength classical chimera 
states undergo transitions from one to multiple domains of incoherence. This additional spatial mod- 
ulation is due to strong coupling interaction and thus cannot be observed in simple phase-oscillator 
models. 
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During recent times the investigation of coupled sys- 
tems has led to joint research efforts bridging between 
diverse fields such as nonlinear dynamics, network sci- 
ence, and statistical physics, with a plethora of applica- 
tions, e.g., in physics, biology, and technology. As the 
numerical resources have developed at fast pace, analysis 
and simulations of large networks with more and more so- 
phisticated coupling schemes have come into reach giving 
rise to an abundance of new dynamical scenarios. Among 
these a very peculiar type of dynamics was first reported 
for the well-known model of phase oscillators. Such a 
network exhibits a hybrid nature combining both coher- 
ent and incoherent parts [H-U, hence the name chimera 
states. The most surprising aspect of this discovery was 
that these states exist in a system of identical oscillators 
coupled in a symmetric ring topology with a symmet- 
ric interaction function. Recent works have shown that 
chimeras are not limited to phase oscillators, but can 
in fact be found in a large variety of different systems. 
These include time-discrete and time-continuous chaotic 
models [H, Q and are not restricted to one spatial di- 
mension. Also two-dimensional configurations allow for 
chimera states [j, 8]- Furthermore, similar scenarios exist 
for time-delayed coupling Q, and their dynamical prop- 
erties and symmetries were subject to theoretical studies 
as well @, [lJJ, [ll| . It was only in the very recent past that 
chimeras were realized in experiments on chemical oscil- 
lators [lH and electro-optical coupled-map lattices [l3| . 
The nonlocality of the coupling - a crucial ingredient for 
chimera states - also suggests an interesting connection 
to material science, see, for instance, magnetic Janus par- 
ticles that undergo a synchronization-induced structural 
transition in a rotating magnetic field 0, HH • Hybrid 



states were also reported in the context of neuroscience 
under the notion of bump states (lfj| . They were later 
confirmed for nonlocally coupled Hodgkin-Huxley mod- 
els fl7l ] and may account for experimental observation of 
par tial synchrony in neural activity during eye movement 

In this letter, we present evidence that the habitat of 
chimeras indeed extends to neural models. This strongly 
suggests their universal appearance. Our findings also 
show that current knowledge of these hybrid states is 
far from being complete: Next to the classical chimera 
state, which exhibits one coherent, phase-locked and one 
incoherent region, we find a new class of dynamics that 
possesses multiple domains of incoherence. 

We consider a ring of N nonlocally coupled FitzHugh- 
Nagumo (FHN) oscillators, whose relevance is not limited 
to neuroscience, but also includes chemical [l9j and opto- 
electronic [20i oscillators and nonlinear electronic circuits 
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where Uk an d Vk are the activator and inhibitor variables, 
respectively [12 , and e > is a small parameter char- 
acterizing a timescale separation, which wc fix at e = 0.05 
throughout the paper. Depending upon the threshold 
parameter ak, each individual FHN unit exhibits either 
oscillatory (\ak\ < 1) or excitable (|ofc| > 1) behavior. In 
this study wc assume that the elements are in the oscil- 
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FIG. 1: (Color online) (a) Limit cycle u = U a (t),v = V a (t) of 
single decoupled FHN oscillator for e — >• 0. (b) Phase interac- 
tion curve H(ip) (red solid) obtained from Eq. © and its ap- 
proximation (blue dashed) by the partial Fourier sum JSJ. (c) 
Phase lag parameter a as a function of control parameters a 
and <j>. The hatched region indicates the parameter range for 
which we find chimera states in simulations of Eqs. {TJ with 
N = 1000, R = 350 and a = 0.05. The thick black curve for 
a = 7r/2 is a guide for the eye. 



latory regime and identical, i.e., au = a G (—1,1). Each 
of them is coupled with its R > nearest neighbors from 
both sides with a coupling strength a > 0. Thus the in- 
dices in Eqs. ((T|) have to be considered modulo N, which 
induces nonlocality in form of a ring topology. This yields 
a coupling radius r — R/N . As we will demonstrate, an 
important feature of Eqs. ([1]) is that it contains not only 
direct u-u and v-v couplings, but also cross-couplings be- 
tween activator (u) and inhibitor (v) variables. For the 
sake of simplicity, we model this feature by a rotational 
coupling matrix 
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depending on a single parameter <j) 6 [— 7r, 7r). Thus we 
can vary four control parameters of different nature: a 
determining the local dynamics, a, i?, and <f> specifying 
the coupling. 

In the following we address the question of appropri- 
ate choice of the coupling phase <j> by a phase-reduction 
technique in the limit e — > 0. In the absence of coupling 
the local dynamics is then periodically oscillating and 
follows a limit cycle (u = U a (t), v = V a (t)), see Fig. HJa), 
with the approximate period [24| 



T = 3 + (l-a 2 )ln 



and two discontinuities at times 
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Here, U a (t) is composed of two distinct solutions of the 
differential equation on the slow activator nullcline [24| 
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The first is defined for t £ [0,ii) and satisfies ini- 
tial condition U a (0) = —2, while the second is defined 
for t G [ti,T) and starts from the point U a (ti) = 2. 
The corresponding slow v- variable is given by V a (t) = 
U a {t) — U^(t)/3. In the weak coupling limit a <C £ <C 1, 
we now perform a phase reduction of Eqs. (JTJ) (25|. Then 
the dynamics of each FHN oscillator is conveniently de- 
scribed by a single scalar phase variables 8k € KmodT, 
which reflects the position of this oscillator along its un- 
perturbed limit cycle. The interaction of these phases, 
to the leading order, is given by a reduced system 
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where the function 
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is T-pcriodic and we use the abbreviations 

p(t) = U a (t) cos <f> + V a (t) sin 4>, 
q(t) = -U a (t) sin</) + V a (t) cos</>, 
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for notational convenience. The function H(tp) typi- 
cally has several discontinuities (see solid (red) curve in 
Fig.fTJb)), because the FHN model becomes a relaxation 
oscillator for e — > 0. 

In order to find suitable parameter values for the gen- 
eration of chimera states, we employ results of Ref. |26| . 
There it was shown that chimera states can be gener- 
ically found in systems of the form of Eq. ([5]) with 
H(tp) = sm(ip + a) if the phase lag parameter a is close 
to but less than 7r/2. Phase interaction curves H(ip) cor- 
responding to a FHN oscillator can be qualitatively ap- 
proximated by their Fourier series truncated at the first 
order 
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FIG. 2: (Color online) (a) Snapshot of the variables Uk for t = 
5000, (b) snapshot in the i)fc)-plane for t = 5000 (black 
lines denote the nullclines of the FHN system), (c) mean phase 
velocities LJk, (d) local order parameter Z^. Parameters: N = 
1000, r = 0.35, cr = 0.1, a = 0.5, <f> = tt/2 - 0.1. 



where the Fourier coefficients ho, h\, and h\ are given by 
h o = ^J o ff(V0#> (9a) 
h i = ^l ffWcos^v)#, (9b) 
h\ = ~J ffty)sinf|^<ty. (9c) 

This yields the following approximate equations for the 
phase-lag parameter a 



^{h\f + (h\f \l(h\f + {h\f 

(10) 

that can be used to pinpoint a region in the parame- 
ter plane (a, </>), which favors the appearance of chimera 
states. Roughly speaking, such states are expected for a 
pronounced off-diagonal coupling (0 « 7r /2), but not for 
a diagonal one [<f> « or </> « 7f) . This prediction is con- 
firmed by numerical simulations, sec the hatched area in 
Fig.QJc). Indeed, if we choose a = 0.5 and cj) = tt/2— 0.1, 
then for cr small enough we obtain a chimera-like solu- 
tion shown in Fig. [2] where we used initial conditions 
randomly distributed on the circle u 2 + v 2 = 4. 

Figure Eta) shows a snapshot of variables Uk at time 
t = 5000. One can clearly distinguish coherent and inco- 
herent parts, a characteristic signature of chimera states. 
Elements that belong to the incoherent part are scattered 
along the limit cycle, see Fig. HJb). The subsystems of 
this region perform a nonuniform rotational motion, but 
neighboring oscillators are not phase-locked. To illustrate 
this, Fig. EJc) shows mean phase velocities for each oscil- 
lator calculated as u>k = ^ttM^/ AT, k = 1, ...,N, where 
Mk is the number of complete rotations around the origin 
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FIG. 3: (Color online) Stability regimes for chimera states 
with one (gray, red), two (dark-gray, green), and three (blue, 
black) incoherent domains in the plane of coupling radius 
r and coupling strength a. Other parameters as in Fig [21 
Hatched regions denote multistable regimes. Insets show typ- 
ical profiles of the mean phase velocities. Black squares (A) 
and circles (B) denote parameter values for the transition sce- 
narios shown in Fig. [4] 



performed by the fc-th unit during the time interval AT. 
The values of uj^ lie on a continuous curve and the in- 
terval of constant u>k corresponds to the coherent region, 
where neighboring elements are phase-locked. This phase 
velocity profile is a clear indication of chimera states and 
similar to the case of coupled phase oscillators. 

The spatial coherence and incoherence of the chimera 
state can be characterized by a real-valued local order 
parameter (El ITT| 
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where Qj = axctaxi(vj /v,j) denotes the geometric phase of 
j-th FHN unit. We use a spatial average with a window 
size of 5 = 25 elements. A local order parameter Zk = 1 
indicates that the fc-th unit belongs to the coherent part 
of the chimera state, and Zk is less than 1 for incoherent 
parts. Figure [21(d) depicts the local order parameter in 
the time interval t 6 [1000,5000], where bright (yellow) 
color denotes the coherent regions. 

For further investigation we fix the values of parame- 
ters a = 0.5 and <fi = tt/2 — 0.1 and vary radius r and 
strength a of the coupling. In previous works, where non- 
locally coupled phase oscillators, discrete maps, Rossler 
and Lorenz systems were considered, chimera states were 
reported for intermediate coupling radii and small cou- 
pling strengths 0, . Figure [3] displays the stability dia- 
gram for chimera states in nonlocally coupled FHN sys- 
tems in a similar parameter range. The gray (red) region 
corresponds to the classical chimera state with one in- 
coherent domain. Surprisingly, for increasing coupling 
strength, we observe qualitatively new types of chimera 
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FIG. 4: (Color online) Transition from a classical chimera state with one incoherent domain to multi-chimera states with 
two (a)-(d), and three (e)-(h) incoherent domains. In each panel the left column shows snapshot of variables Uk, and the right 
column shows the corresponding mean phase velocities. Coupling radius and strength follow the black squares (A) in Fig.[3]for 
panels (a)-(d) and black circles (B, a — — 3r + 1) in Fig. [3]for panels (e)-(h). Other parameters as in Fig. [2] 



states, which have two or even three incoherent domains 
indicated by the dark-gray (green) and black (blue) re- 
gions, respectively. We call these states multi-chimera 
states. Their additional spatial modulation cannot be ex- 
plained in terms of phase interaction only Q. Near the 
borders of the different regimes, multistability is found 
and indicated by hatching. In this region of coexistence 
of one and two (or one and three) incoherent domains 
the particular realization depends on the choice of initial 
conditions. 

The transition from the classical to the multi-chimera 
state is illustrated in Fig. 01 Following the lines with 
black squares (A) and black circles (B) in Fig. H 
Figs. SJa)-(d) and S£e)-(h) demonstrate how multi- 
chimeras with two and three incoherent domains develop 
from the classical chimera with one incoherent part. The 
velocity profile shows a dip, which becomes more pro- 
nounced as a increases and eventually reaches down to 
the level of the coherent part. For smaller coupling 
radii, Eq. ([1]) exhibits increasing multistability. There 
wc also find chimera states with more than three incoher- 



ent parts, but these states have relatively small stability 
domains and are visible for large system sizes N only. 

In conclusion, we have reported the existence of 
chimera states for relaxation oscillators, which are of 
slow-fast type and exhibit a timescale separation between 
activator and inhibitor. Applying a phase-reduction tech- 
nique, wc have identified nonlocal off-diagonal coupling 
to be a crucial ingredient for the occurrence of these 
hybrid states that exhibit coexistent coherent and in- 
coherent domains. Furthermore, we have found a new 
type of multi-chimera states that consist of multiple do- 
mains of incoherence. They appear as a result of strong 
coupling interaction and thus cannot be found in simple 
phase models. Our findings corroborate the universal ex- 
istence of chimera states, which were previously reported 
for phase oscillators and time-discrete or time-continuous 
chaotic systems. 
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